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Antisense RNA (asRNA) COOLAIR is expressed at A. thaliana FLOWERING LOCUS C (FLO in 
response to winter temperatures. Its contribution to cold-induced silencing of FLC was 
proposed but its functional and evolutionary significance remain unclear. Here we identify a 
highly conserved block containing the COOLAIR first exon and core promoter at the 3' end of 
several FLC orthologues. Furthermore, asRNAs related to COOLAIR are expressed at FLC loci 
in the perennials A. alpina and A. lyrata, although some splicing variants differ from A. thaliana. 
Study of the A. alpina orthologue, PERPETUAL FLOWERING 1 {PEPl), demonstrates that 
AaCOOLAIR is induced each winter of the perennial life cycle. Introduction of PEPl into 
A. thaliana reveals that AaCOOLAIR c/s-elements confer cold-inducibility in this heterologous 
species while the difference between PEPl and FLC mRNA patterns depends on both 
c/s-elements and species-specific frans-acting factors. Thus, expression of COOLAIR is highly 
conserved, supporting its importance in FLC regulation. 
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Many plant species growing in temperate climates flower 
in response to prolonged low winter temperatures, 
a process called vernalization^'^. In Arabidopsis 
thaliana FLOWERING LOCUS C (PLC) plays a central role in 
vernalization requirement, encoding a MADS box transcription 
factor that delays flowering until the plant is exposed to low 
temperatures^'"*. On extended exposure to cold PLC mRNA 
level is progressively reduced, so that after several weeks it 
reaches trough levels^. On return to warmer temperatures, 
mimicking spring conditions, PLC repression is maintained and 
flowering occurs. 

Repression of PLC by cold involves epigenetic regulation. 
Histone modifications, particularly trimethylation of lysine 27 on 
histone 3 (H3K27me3), increase at PLC during vernalization 
through the activity of Polycomb Repressive Complex 2 (PRC2) 
and later spread across the gene^'^. This mechanism requires 
accessory components of the PRC2 that are cold-induced, such as 
PHD-finger protein VERNALIZATION INSENSITIVE 3 
(VIN3)^"^. VIN3 as well as core PRC2 proteins such as 
VERNALIZATION 2 (VRN2) are required to maintain PLC in 
a repressed state after vernalization^'^ '^^. 

Repression of PLC is also proposed to involve long non- coding 
RNAs (IncRNAs) expressed at the locus in response to cold*^ A 
sense IncRNA called COLD AIR is expressed during vernalization 
from a region in the long first intron*^ termed the vernalization 
response element (VRE)*^. The VRE is downstream of the 
nucleation point at which chromatin modifications initially 
occur^ and was previously shown to be important for stable 
silencing of PLC COLD AIR binds PRC2 components and 
promotes PLC stable repression*^. Earlier during vernalization, 
before expression of COLD AIR, a set of alternatively spliced and 
polyadenylated antisense RNAs (asRNAs) termed COOLAIR are 
transiently expressed from the 3' end of FLC*^~*^. Fusion of the 
COOLAIR promoter at the 3' end of a green fluorescent protein 
(GFP) transgene can cause partial repression of GFP expression 
during vernalization*^. COOLAIR has been shown to participate 
in PLC downregulation by the autonomous pathway but its 
function in the cold-induced silencing of PLC is still 
unclear *^'*^'*^. Indeed, PLC marker genes lacking the 
COOLAIR promoter silence in cold*"* and PLC insertion alleles 
from which COOLAIR expression is not detectable are also 
silenced by vernalization*^. Recently, COOLAIR transcription 
was shown to be downregulated by the formation of a RNA- 
DNA hybrid within its promoter. The so-called R-loop is 
stabilized by the homeodomain protein AtNDX (ref. 20). 

PLC homologues are found throughout the Brassicaceae family, 
of which A. thaliana is a member. Expression of some of these 
homologues has been shown to be reduced by vernalization^ *"^^ 
and genetic evidence demonstrates that some act as floral 
repressors^ *'^^. In addition to their role in vernalization, 
differences in expression of PLC orthologues can contribute to 
diversification of life history within the Brassicaceae. A. thaliana 
is a representative monocarpic annual species of the family that 
flowers only once and then dies after seed set. However, other 
members of the family are polycarpic perennials and undergo 
multiple reproductive episodes during their life. Interestingly, in 
the perennial species Arabis alpina, A. halleri, Cardamine flexuosa 
and A. lyrata the expression patterns of PLC orthologues differ 
from A. thaliana PLC. The PLC genes in the perennials are not 
stably repressed by vernalization but their expression increases 
after return to warm temperatures^ *'^^'^^'^^. The A. alpina PLC 
orthologue, PERPETUAL PLOWERING 1 (PEPl), has additional 
perennial specific functions. For example, the increase in PEPl 
mRNA after vernalization confers the return to vegetative growth 
after flowering until plants are exposed to low temperatures the 
following winter^*. 



Comparison of transcriptional regulatory regions from several 
closely related species has been employed to identif)^ conserved 
c/s-regulatory elements in different taxa^^~^* including 
plants^^"^"*. In contrast, the extent to which IncRNAs diverge 
during evolution has not been extensively analysed in plants, 
although in other systems such as mammals they appear to evolve 
rapidly between species^^'^^. 

To better understand what drives the conserved cold- regula- 
tion of PLC orthologues and to assess the functional significance 
of IncRNA at this locus, we combine phylogenetic shadowing of 
PLC c/s-regulatory sequences with structural and temporal 
analysis of antisense transcription, especially at PEPU the best- 
characterized PLC orthologue. Our data reveal that the 3' end 
sequences at PLC orthologues are more conserved than the 5' 
promoter regions and that they promote antisense transcription 
at A. alpina PEPl and A. lyrata PLC genes similar to what 
was observed in A. thaliana PLC. Introduction of PEPl into 
A. thaliana shows that it contains all the c/s-acting sequences 
required for cold repression and COOLAIR expression in this 
heterologous species. In contrast, the reactivation of PEPl after 
vernalization seems to involve different trans -acting factors 
between A. alpina and A. thaliana. Evolutionary conservation 
of COOLAIR and its seasonal expression in perennial A. alpina 
supports a significant role for this asRNA in repression of PLC 
gene expression during vernalization. 

Results 

Conserved sequence blocks at the 3' end of PLC orthologues. In 

A. thaliana cold-induced repression of PLC mRNA requires c/s- 
regulatory sequences in the promoter and first intron*^'*"*. 
Whether these c/s-regulatory elements are evolutionarily 
conserved was tested by comparing genomic sequences of PLC 
orthologues from a range of Brassicaceae species (Methods; Fig. 1; 
Supplementary Figs l-3).With the exception of PEPl in A. 
alpina, PLC orthologues from A. lyrata, Capsella rubella and 
Thellungiella halophila exhibit the same gene structure as A. 
thaliana PLC. A. alpina PEPl harbours a tandem duplication 
(PEP la and PEP lb) containing part of the promoter, the first 
exon and part of intron 1 (Fig. 2a; ref. 21). The two copies of exon 
1 have individual transcriptional start sites and each is spliced to 
the unique exon 2 giving rise to two overlapping transcripts^^. 
Two tandemly arranged full-length copies of PLC, called AlPLCl 
and AIPLC2, are present in A. lyrata, and both are 
transcriptionally active^^'^^. Promoter sequences between the 
translational start site (ATG) of PLC and the end of the next 
upstream gene were compared revealing two main blocks of 
homology (Fig. la and Supplementary Fig. 1). The first block was 
~ 400 bp long and included the proximal promoter as well as the 
5' untranslated leader. In A. alpina, homology in this region was 
reduced and mainly restricted to the untranslated leaders of 
PEP la and PEP lb (Supplementary Fig. Ig). The second region 
was ~ 2 kb further upstream, and in A. alpina was found only in 
PEP la. Apart from these two conserved blocks, the PLC promoter 
has diverged rapidly during evolution, even between the closely 
related species A. thaliana and A. lyrata, as observed previously^ . 

Alignment of the intron 1 sequence of PLC with those of the 
orthologues showed that the extensive homology present in the 
closest species decreases in the more distantly related species until 
it is restricted to each end of the first intron and to the VRE 
(Supplementary Fig. 2a-e). These regions are important for PLC 
regulation by cold*^'*"*. Sequences of the PLC orthologues 
showing homology to A. thaliana VRE were referred to as 
VRE-like and were found in all species examined (Fig. lb). 
Interestingly, A. alpina PEPl carries two VRE-like copies in 
intron 1, AaVRE-like 1 and AaVRE-like 2, both of which are 
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Figure 1 | Conservation and divergence of non-coding regulatory c/s-elements annong F/.C orthologues. (a) mVista pair-wise alignments of the 5' intergenic 
regions of FLC orthologues. From top to bottom, A. lyrata FLO, A. lyrata FLC2, C. rubella FLC, T. halophila FLC and A. alpina PEPIa aligned to A. thaliana FLC. 
Alignment of PEPIb and A. thaliana FLC is presented in Supplementary Fig. If. (b) GATA pair-wise alignment of A. thaliana VRE to the first intron of each of 
the FLC orthologues mentioned in a. Regions showing homology to the VRE are annotated as VRE-like. For A. alpina PEPl, VRE-like sequences were detected 
only in the region spanning exon lb and exon 2 (intron 1 of PEPIb). Scale bar, 1 kb. (c) mVista pair-wise alignment of the 3' regions of FLC orthologues. PEPl 
3' end was aligned to the 3' intergenic regions of A. thaliana FLC, A. lyrata FLCl, A. lyrata FLC2, C rubella and J. halophila FLC orthologues. Coloured areas in 
mVista graphs illustrate stretches of homology greater than 75% identity at the nucleotide level. Pink, regions of homology in the promoter region (a) or at 
the 3' end (c); dark blue, exonic sequences; light blue, untranslated region. Grey and pink boxes in GATA plots are homologous regions and inverted 
homologous regions respectively (the darker the box, the higher the similarity). 



located downstream of exon lb. The COLD AIR promoter is 
contained within the most S' 100 bp of the VRE in 
A. thaliana^'^'^^ , and was the least conserved part of the VRE- 
like sequences, showing no detectable conservation even between 
A. thaliana and A. lyrata (Fig. lb and Supplementary Fig. 2a,b). 

The FLC 3' end harbours a promoter that drives the expression 
of the non- coding COOLAIR asRNA^^ Alignment of the ?>' end 
of FLC to those of the orthologues showed strong conservation of 
a 150 bp block after the polyadenylation site (Fig. Ic and 
Supplementary Fig. 3). Strikingly, this ?>' region is more widely 
conserved among the analysed species than either the S' 
intergenic region or first intron sequences implicated in FLC 
transcriptional regulation. In A. thaliana this evolutionarily 
conserved block is present within the COOLAIR first exon, the 
region bound by AtNDX to stabilize the R-loop and regulate 
COOLAIR expression^^. Conservation of the ?>' ends of 
orthologous plant genes has not been extensively analysed, so 
to determine whether the conservation at this location of FLC is 
unusual, 10 other genes encoding MADS box transcription 
factors implicated in reproductive development were examined. 
The intergenic regions at the ?>' ends of orthologues from 
A. lyrata, C rubella, T. halophila, A. alpina and A. thaliana were 
aligned (Supplementary Fig. 4). For nine of the genes studied, 
conservation was observed only in the transcribed untranslated 
region present at the ?>' end of the cDNA, and this reduced steeply 
outside of the transcribed region, especially in more distantly 
related species. Therefore, strong conservation after the 3' end of 



the transcribed region, similar to that present in FLC, is not 
detected among these genes. However, the remaining gene, 
SHATTERPROOF! {SHP2), exhibited conservation beyond the 3' 
UTR that was similar to that observed for FLC. Interestingly, a 
potential natural antisense IncRNA is also annotated at its 3' end 
and named At2g42835 (Supplementary Fig. 4). 

Antisense transcription occurs at PEPl 3' end. The highly 
conserved 150-bp block present at the 3' end of all FLC ortho- 
logues suggested that asRNAs might also be transcribed at these 
genes. The presence of antisense transcription was therefore 
tested at A. alpina PEPl, the most characterized FLC orthologue. 
Several forms of PEPl asRNA were detected by y tiling' RT-PCR 
and 3' RACE PCR (see Methods and Fig. 2a). The detected 
asRNAs were assigned to classes according to their structures and 
named based on the nomenclature used for A. thaliana COOL- 
AIR^ ^. Class I and II A. alpina PEPl asRNAs have similar 
structures to those of A. thaliana COOLAIR^^'^^'^^, and were 
named AaCOOLAIR I and AaCOOLAIR II, respectively. Class I 
asRNAs are polyadenylated within intron 6 of PEPl, while those 
in Class II are spliced to a distal location in the PEPl promoter. 
We further tested whether AaCOOLAIR II are polyadenylated 
within the promoter of PEP la or PEP lb by making use of the 
accession Dorfertal (Dor), which contains polymorphisms in 
these regions^^. Performing '3' tiling' RT-PCR for Dor PEPl 
showed that class II asRNAs terminate in the promoter regions of 
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Figure 2 | Cold-induced antisense non-coding RNAs are transcribed at the PEP1 locus in A. alpina and follow a seasonal pattern, (a) Class I (AaCOOLAIR I), 
class II (AaCOOLAIR II) and class III (AaCOOLAIR III) antisense RNAs transcribed at PEP1 locus. Black boxes are exons; lines are introns and other 
non-coding regions. Grey boxes indicate the region that is duplicated in tandem. Arrows indicate transcriptional start sites at PEPT. Dots show PEP1 
antisense transcript polyadenylation sites; blue, Class I; light orange, Class III; dark orange. Class II. Asterisks indicate the most abundant form among 
splicing variants obtained by 3' RACE (Classes I and III) or RT-qPCR product sequencing (Class II). Scale bar, 1 kb. (b) AaCOOLAIR I, II and III 
cold-inducibility. RT-qPCR on leaf RNA of A alpina plants exposed to 0, 7, 21, 35 and 84 days of vernalization (VO, V7, V21, V35 and V84) followed by 7, 21 
and 35 days of growth in normal growth temperatures (V + 7, V + 21 and V + 35). Transcript levels were expressed relative to those of the reference gene 
RAN3±s.di. (n = 3 technical replicates), (c) AaCOOLAIR I, AaCOOLAIR III, PEPl and AaVIN3 expression patterns in A. alpina leaves during the vernalization 
time course shown in b. Transcripts levels were measured by RT-qPCR relative to those of the reference gene RAN3±s.6. (n = 3 technical replicates). 
The maximum level of expression of each gene is set at 100%. (d) Seasonal expression patterns of PEPl, AaCOOLAIR I, AaCOOLAIR III and AaVINS 
over two successive vernalization treatments. A. alpina plants were vernalized for 0, 7, 21, 35 and 84 days (VO, V7, V21, V35 and V84) followed by 7, 
21 and 35 days of growth in normal temperature (V + 7, V + 21 and V + 35), then vernalized again for 7, 21, 35, 56 and 84 days (VV7, VV21, VV35, 
VV56 and VV84) and grown for 7, 21 and 35 more days in normal growth temperatures (VV + 7, VV + 21 and VV + 35). Grey areas indicate 
the cold treatments. Transcript levels were measured by RT-qPCR relative to those of the reference gene RAN3±s.6. (n = 3 technical replicates). 
The maximum level of expression of each gene is set at 100%. A biological replicate is shown in Supplementary Fig. 8. 



both PEP la and PEP lb (Supplementary Fig. 5). A third class of 
PEPl asRNA (AaCOOLAIR III) not documented in A. thaliana 
was also isolated, and is polyadenylated in intron 1 in a region 
that shows no homology with PLC (Fig. 2a and Supplementary 
Fig. 2e). We further confirmed the conservation of antisense 
transcription at the ?>' end of PLC orthologues by amplif)^ing 
asRNA generated at A. lyrata FECI and FECI genes, both of 
which show homology to COOLAIR (Fig. Ic). Class I but not 
class II asRNA could be detected and harbour the typical splicing 
events of COOLAIR I asRNAs (Supplementary Fig. 6a). 

In contrast to the detection of asRNAs, non-coding sense RNA 
corresponding to COLD AIR could not be amplified at PEPl 
after 5' RACE PCR using several sets of primers downstream of 
the VRE-like 1 and the VRE-like 2 in intron 1 {5' RACE 'Tiling 
PCR' in Supplementary Methods, Supplementary Table 1). 
Vernalized A. thaliana FRI EEC plants were also used as a 



control but the previously described COLD AIR transcript could 
not be reliably detected, suggesting that our vernalization 
conditions or methodology may not be optimal for the 
detection of these non-polyadenylated transcripts. Therefore, 
COLD AIR could not be detected in A. alpina, but its presence 
cannot be excluded. 



Induction of AaCOOLAIR by cold in A, alpina. The levels of 
the different AaCOOLAIR forms and their cold-inducibility were 
tested by RT-qPCR in A. alpina leaves during vernalization 
(Fig. 2b). AaCOOLAIR I and III were the most abundant forms. 
They showed two peaks of expression during cold exposure, after 
21 and 84 days in the cold, and fell drastically when the plants 
were returned to higher temperatures after vernalization. Cold- 
inducibility seems to be reminiscent of class I asRNA at EEC 
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orthologues as it was also observed for A. lyrata AlCOOLAIR I 
(Supplementary Fig. 6b). By contrast, class II asRNA could not be 
reliably amplified in A. lyrata. AaCOOLAIR II was detected but 
expressed at a very low level with no significant change in 
abundance during vernalization (Fig. 2b). 

In A. thaliana transcription of VIN3 is also induced during 
vernalization and is required for epigenetic silencing of FLC^. The 
A. alpina VIN3 orthologue was isolated based on sequence 
homology and showed conservation of synteny with A. thaliana 
(Supplementary Fig. 7). AaVIN3 mRNA increased during cold, as 
observed in A. thaliana, and fell steeply at the end of vernalization 
(Fig. 2c). 

Seasonal cycling of AaVIN3 mRNA and AaCOOLAIR in A. alpina. 

PEPl expression repeatedly responds to seasonal changes in 
temperature contributing to the perennial life cycle of A. alpina^^. 
Whether AaCOOLAIR I and III as well as AaVIN3 have similar 
responses was tested. A. alpina plants were exposed to two 
successive vernalization treatments separated by 35 days at normal 
growth temperatures. Leaves were harvested at regular time 
intervals during the treatment (see Methods). As expected, PEPl 
mRNA levels responded to changes in temperature, being 
downregulated during cold and upregulated when plants were 
transferred to warm temperatures (Fig. 2d, Supplementary Fig. 8). 
Interestingly, abundance of AaCOOLAIR I and III cycled inversely 
to PEPl mRNA supporting the idea that they play a role in PEPl 
repression. The timing of their maximal expression differed 
between the two cold treatments so that they reached peak level 
earlier in the second treatment than the first. AaVIN3 mRNA 
levels were dynamically induced by cold during both vernalization 
treatments presumably promoting PEPl repression and the 
flowering of the plants during each following warm period 
(Fig. 2c,d, Supplementary Fig. 8). Interestingly, the induction was 
more than 10-fold stronger during the second cold treatment than 
the first. Therefore, expression of AaCOOLAIR and AaVIN3 
mRNA might be primed by the first vernalization treatment to 
respond more strongly during the second. 

PEP lb is stably silenced by cold in transgenic A. thaliana. To 

investigate whether the molecular mechanisms regulating 
expression of FLC orthologues are evolutionarily conserved, A. 
alpina PEPlb (Fig. 3a), the most highly expressed A. alpina PEPl 
gene^^ and the one containing the 3' sequences encoding 
AaCOOLAIR, was transferred into A. thaliana Columbia FRI 
flc-3 (ref. 3). Transgenic FRI flc-3 plants containing the PEPlb 
genomic fragment displayed a small but statistically significant 
(P< 0.005) delay in flowering compared with the FRI flc-3 
progenitor line (Fig. 3b, Supplementary Table 2). By contrast, the 
control line FRI FLC that carries an active FLC gene flowered 
much later than FRI flc-3 (Fig. 3b). PEPl mRNA was also 
expressed in FRfi^ flc-3 PEPlb transgenic lines at lower levels 
than FLC mRNA in the FRI FLC control plants (Fig. 3c, VO). 
These data suggest that the diverged PEPlb promoter might not 
be efficiently recognized in A. thaliana and cannot promote PEPl 
mRNA at the threshold level necessary to strongly delay flowering 
time. To test the effect of cold on PEPl mRNA level transgenic 
plants were vernalized for 40 days and subsequently transferred 
to warm temperatures for 20 days (Fig. 3c). PEPl mRNA levels 
were reduced during vernalization as observed for FLC mRNA in 
FRI FLC control plants (V40, Fig. 3c). Interestingly, in this 
heterologous species PEPl b mRNA level remained low when the 
transgenic plants were returned to normal growth temperatures 
after vernalization, similar to what is observed for FLC stable 
repression (V40+10, V40 + 20; Fig. 3c and Supplementary 
Fig. 9a) and in contrast to the unstable repression of PEPl 
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Figure 3 | Transcriptional activity of PEPlb in A. thaliana FRI flc-3 PEPlb 
transgenic plants, (a) Schematic representation of PEPl locus represented 
as described in Fig. 2a. The rectangle indicates the extent of 13.4 kb PEPlb 
transgene transferred into A. thaliana FRI flc-3 plants. Scale bar, 1 kb. 
(b) Flowering tinnes of A. thaliana FRI flc-3 PEPlb independent trans- 
formants, FRI flc-3 parental line and FRI FLC control line as measured by 
rosette leaf number (RLN) at bolting ± s.d. (n = 9-23 individuals). Asterisks 
indicate flowering time statistically different from the one of FRI flc-3 
plants (Student's f-test, P value < 0.005). (c) PEPlb mRNA levels in 
A. thaliana FRI flc-3 PEPlb independent transformants and FLC mRNA 
levels in the parental FRI ^'^ flc-3 line and the control FRI ^'^ FLC line. Time 
points are before vernalization (VO), after 40 days of vernalization (V40), 
as well as after 40 days of vernalization followed by 10 or 20 days of 
growth at normal growth temperatures (V40 + 10 and V40 + 20). 
Transcript levels were measured by RT-qPCR relative to those of the 
reference gene ACTIN±s.6. (n = 3 technical replicates). Panels b and c 
share colour-code legends. 



found in A. alpina. However, the low level of expression of PEPlb 
in A. thaliana might also reduce the sensitivity of detection of 
reactivation after vernalization. Despite this reservation, the 
species-specific variation in trans-acting regulation is likely 
responsible for the differences in the stability of PEPl 
repression observed between A. thaliana and A. alpina. 

Epigenetic regulation of PEPl transgene in A. thaliana. Whe- 
ther A. alpina AaCOOLAIR was generated by mechanisms con- 
served in A. thaliana was examined. AaCOOLAIR was detected 
in the FRI flc-3 PEPlb transformants during vernalization 
(Fig. 4a), suggesting that trans-acting 
thaliana can promote AaCOOLAIR 



factors present in A. 
expression, including 
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Figure 4 | Epigenetic control of PEPIb in Arabidopsis FRI ^'^ flc-3 PEPlb 
transgenic plants, (a) PEPlb, AtVIN3 and AaCOOLAIR I and III expression in FRI 
PEPlb 18-3 transformant before vernalization (VO) after 10, 20, 30 or 
40 days of vernalization (V10, V20, V30 or V40) and after 40 days of 
vernalization followed by 10 or 20 days of growth in normal temperature 
(V40 + 10, V40 + 20). (b) FLC, AtVIN3, AtCOOLAIR I and II expression in FRI 
FLC control plants during the vernalization time course described in a. In a 
and b mRNA levels were measured by RT-qPCR relative to those of the 
reference gene ACTINls.d. (n = 3 technical replicates) and presented relative 
to the maximum level of expression of each gene, which is set as 100%. 
(c) PEPlb mRNA level before vernalization (VO) and after 40 days 
vernalization followed by 20 days of growth at normal temperature 
(V40 + 20) in FRI flc-3 PEPlb 18-3 parental line and FRI flc-3 PEPlb 18-3 
vrn2-l fca-1 plants. PEPlb mRNA level is compared with FLC mRNA in FRI FLC, 
FRI flc-3 and wn2-l fca-1 plants. mRNA levels were measured by RT-qPCR 
relative to those of the reference gene ACTIN ± s.d. (n = 3 technical replicates). 

AaCOOLAIR III that is not found at FLC. In the control Une FRI 
FLQ AtCOOLAIR I and AtCOOLAIR II were induced after 20 
days in cold followed by accumulation of AtVIN3 mRNA and 



stable repression of FLC transcripts (Fig. 4b). Similarly, in 
transgenic FRI flc-3 PEPlb plants AaCOOLAIR I and 
AaCOOLAIR III were also induced by cold and their patterns of 
expression resembled those of A. thaliana COOLAIR I and II 
(Fig. 4a,b). Therefore, AaCOOLAIR transcripts are probably 
controlled by c/s- regulatory elements present in PEPl that are 
recognized by trans-dLCting factors conserved between A. thaliana 
and A. alpina. 

In A. thaliana^ FLC is stably silenced during cold by the 
deposition of chromatin- repressive marks (H3K27me3) and their 
spreading over the locus after return to normal growth 
temperature^'^^"'*^ VRN2, a component of PRC2, is required 
for both establishing and maintaining this epigenetic silenced 
state ^^'"^^ To test whether chromatin modification by PRC2 is 
responsible for the stable repression of the PEPlb transgene 
observed in Arabidopsis FRI flc-3 PEPlb we investigated the 
effect of a vrn2 mutant allele on PEPlb expression. The vrn2-l 
fca-1 mutations were introgressed into FRI flc-3 PEPlb 
background (see Methods). The fca-1 mutation impairs the 
autonomous flowering pathway and therefore increases FLC 
expression before vernalization but does not affect stable 
repression of FLC after vernalization FLC and PEPl mRNA 
levels were measured in control and introgressed plants before 
and after 40 days of cold treatment (Fig. 4c, Supplementary 
Fig. 9a). The vrn2-l fca-1 mutants failed to stably repress FLC 
mRNA expression upon vernalization, showing much higher FLC 
levels than in the FRf'^ FLC control line 20 days after return to 
warm temperature (Fig. 4c, V40 + 20). Such deregulation could 
also be observed for the residual endogenous FLC expression in 
FRI flc-3 PEPlb vrn2-l fca-1 (Supplementary Fig. 9b, 
V40 + 20). PEPl mRNA levels in FRI flc-3 PEPlb vrn2-l fca- 
1 transgenic lines were higher before vernalization compared with 
FRI flc-3 PEPlb 18-3 parental line (Fig. 4c, VO), similar to what 
was described for the endogenous FLC in vrn2-l fca-1 mutant^^. 
Moreover, PEPl showed a tremendous increase of expression 20 
days after return to warm temperature in this background, 
whereas its mRNA levels were almost undetectable in FRI flc-3 
PEPlb 18-3 parental line in which VRN2 is active (Fig. 4c, 
V40 + 20). These data suggest that the stable repression of PEPlb 
observed in A. thaliana FRI flc-3 PEPlb transgenic lines is due 
to epigenetic silencing promoted by PRC2 either directly on 
PEPlb or indirectly by affecting expression of a trans-acting 
regulator. 



Discussion 

In this study, we demonstrate strong conservation of asRNAs 
expressed at FLC orthologues during vernalization in the related 
species A. alpina and A. lyrata. In perennial A. alpina we find that 
these asRNAs cycle in response to temperature changes through 
repeated vernalization treatments in a strictly inverse relationship 
to PEPl mRNA, which argues for their involvement in the 
repression of FLC genes by cold. Phylogenetic shadowing 
combined with interspecies gene transfer showed that sequences 
required for transcriptional activation of FLC orthologues in 
A. thaliana and A. alpina have probably diverged, while variation 
in the stability of repression of these genes by vernalization is 
probably caused by trans-acting differences. These results provide 
insights into evolution of the structure and regulation of FLC, a 
gene that controls intra- and inter- species variation in flowering 
behaviour in Brassicaceae species. 

In A. thaliana and A. alpina, the orthologous genes FLC and 
PEPl show different responses to cold. FLC is stably repressed by 
vernalization, whereas cold repression of PEPl is unstable so that 
PEPl mRNA levels increase again after return to warm 
conditions^ ^ This difference in the stability of repression was 
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not observed in transgenic A. thaliana carrying a genomic 
fragment of PEP lb, in which PEP lb mRNA stayed low after 
return to warm temperature. Introduction of a vrn2 mutation 
released PEP lb repression after return to warm conditions, 
suggesting that PEP lb stable repression is not caused by an 
inability to reactivate transcription. Thus, PEP lb must contain all 
the sequences required for an active process of stable repression 
in A. thaliana. The VRE and the COLDAIR IncRNA transcribed 
from it were previously shown to be required for stable silencing 
of PLC in A. thaliana^^-^^. We have shown that the VRE-like 
sequences present in PEP lb lack homology to the COLDAIR 
promoter. This divergence however does not impair PEP lb stable 
silencing in A. thaliana by a mechanism involving the PRC2 
component VRN2. However, we cannot exclude that the effect of 
PRC2 on PEP lb in A. thaliana is indirect as H3K27me3 levels at 
the gene were not measured. The unstable PEPl repression 
observed in A. alpina is most likely due to differences in the 
activity of trans- acting factors compared with A. thaliana and not 
to c/s-acting differences such as those in the VRE, because PEPlb 
was stably repressed in A. thaliana, although the low level of 
expression of PEPlb in A. thaliana might also compromise our 
ability to detect reactivation. Trans-acting factors required for 
stable repression of PLC such as the components and accessory 
proteins of PRC2 (refs 7,8,10,41), LIKE HETEROCHROMATIN 
PROTEIN 1 (LHPl) or VERNALIZATION 1 (refs 13,42,43) are 
present in the A. alpina genome. Their absence, therefore, does 
not explain the instability of the H3K27me3 repressive marks at 
PEPl (ref. 21). An alternative possibility is that A. alpina encodes 
enzymes that actively remove H3K27me3 from PEPl allowing its 
reactivation after vernalization and that this activity is absent 
from annual A. thaliana. H3K27me3 demethylases are present in 
A. thaliana"^"^, and members of this enzyme family might 
specifically target PEPl after vernalization in A. alpina. 

PEPlb was only expressed at low levels in non-vernalized 
transgenic A. thaliana, suggesting that c/s-regulatory elements 
required for PEPlb expression prior vernalization in A. alpina are 
not efficiently recognized in the heterologous species. Indeed, we 
found little 5' homology between PLC and PEPl, and what is 
present is restricted to a proximal block in PEP la and PEPlb as 
well as a distal block only upstream of PEP la. The latter is located 
2kb upstream of PLC translational start site in A. thaliana and 
outside the region shown to be sufficient for PLC regulation 
Other promoter regions moreover have been implicated in 
activation and repression of PLC expression but did not exhibit 
high conservation in the PEPl promoter. RNA or RNA 
polymerase-associated complexes have been described to be 
positive regulators of PLC expression"*^""*^. These appear to act 
close to the transcriptional start site, which may explain the 
homology detected around this region. The proposed binding site 
of one of them, SUF4 (ref. 47), is however not conserved in PEPl. 
Also, one of two populations of A. lyrata studied by Kemi et al.^^ 
was monomorphic for a deletion in FLCl promoter containing 
this binding site of SUF4. These findings further suggest 
that different mechanisms might promote PEPl expression in 
A. alpina as well as FLCl expression in A. lyrata. 

Our analysis demonstrated that the first introns of PLC 
orthologues, especially the flanks and the VRE-like sequences that 
are involved respectively in activating expression of PLC and 
stably repressing the gene in response to vernalization were 
generally more conserved than the promoter regions. In PEPl 
these sequences might also contribute to its downregulation by 
cold although the VRE-like sequences present in the first intron 
do not confer stable repression in A. alpina. 

Both sense and antisense non- coding RNAs have been reported 
at FLC^^'^^ Similar to the PLC non-coding sense RNA 
COLDAIR, the HOTAIR non-coding sense RNA of humans 



was proposed to recruit PRC2 to the HOX gene cluster"^^'"^^. 
However, sequence comparisons and deletion of the mouse 
HOTAIR gene indicate that this regulatory function has diverged 
between humans and mouse^^. Our phylogenetic shadowing 
approach showed that the COLDAIR promoter described in 
A. thaliana is not strongly conserved even in annual species in 
which PLC orthologues are stably repressed by cold. Additionally, 
we could not detect the sense COLDAIR RNA either in A. alpina 
or A. thaliana. These data raise the question whether this sense 
IncRNA is conserved in the Brassicaceae, but the absence of a 
reliable positive control prevents us from drawing a definitive 
conclusion. 

By contrast cold- inducible asRNAs were readily detected at 
PLC genes of A. lyrata and at PEPl of A. alpina, respectively the 
closest and the most distantly related species to A. thaliana that 
we examined. The COOLAIR proximal promoter and the 
COOLAIR first exon at the 3' end of all PLC orthologues were 
more highly conserved than most of the promoter region or first 
intron regulatory sequences of PLC. Such a high conservation was 
not observed at the 3' ends of other MADS box genes, with the 
exception of SHP2 for which a potential natural antisense 
IncRNA is also annotated after the 3' UTR. These analyses suggest 
that COOLAIR IncRNAs are under high constraint and likely 
have a role in cold repression of PLC. However, PLC derivatives 
lacking the 3' end region or in which insertions prevent 
expression of COOLAIR still silence in the cold^"*'^^, suggesting 
that COOLAIR is not essential for PLC silencing. Functional 
redundancy between COOLAIR and the VRE/COLDAIR might 
prevent detection of an effect of COOLAIR. Deletion of the VRE/ 
COLDAIR prevents silencing but does not affect cold-mediated 
repression of PLC mRNA level^^, so it would be interesting to test 
whether loss of both VRE/COLDAIR and COOLAIR prevents 
PLC downregulation in the cold. Nevertheless, current studies 
indicate that COOLAIR is not absolutely required for PLC cold- 
mediated repression or silencing, but might rather quantitatively 
enhance one or both of these processes. COOLAIR might also 
contribute more to PLC regulation in natural conditions where 
temperature and other environmental factors fluctuate widely. 
Mutations affecting flowering time were indeed previously shown 
to have different effects in natural conditions than in growth 
chamber experiments^^. 

Studies of the A. thaliana fpa mutant, which is affected in 
polyadenylation signal recognition, suggest that the distal form of 
asRNA (COOLAIR II) might not play a role in PLC repres- 
^[qyi^6,17 q^j. observation that the proximal COOLAIR I is the 
most conserved asRNA between A. thaliana, A. lyrata and 
A. alpina would support this hypothesis. The presence of distal 
asRNAs (AtCOOLAIR II, AaCOOLAIR II and III) may only be a 
consequence of inefficient polyadenylation at the termination site 
for Class I COOLAIR. Taken together our data suggest that Class 
I COOLAIR asRNAs are those with the major biological function 
and that their increased abundance during vernalization in 
A. thaliana, A. lyrata and A. alpina is a widely conserved feature 
associated with repression of PLC genes in cold. 

The polycarpic, perennial life cycle of A. alpina causes the plant 
to be exposed to vernalization many times during its life time, 
whereas annual, monocarpic A. thaliana is exposed to vernaliza- 
tion only once. We observed that PEPl mRNA cycled with 
temperature, so that it was reactivated in the warm after 
vernalization^ ^ and repressed again by a second vernalization 
treatment. AaCOOLAIR and AaVIN3 cycled in an inverse pattern 
to that of PEPl. Interestingly, they appeared to be primed by the 
first vernalization treatment, resulting in a more rapid and 
enhanced rise in their expression during the second vernalization 
treatment. These differences suggest that perennial plants can 
maintain a memory of their first winter. Molecular mechanisms 
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by which this occurs are unknown but might involve epigenetic 
changes or changes in chromosomal configuration. Such 
mechanisms are unlikely to involve increased H3K27me3 levels 
at PEPl because this chromatin mark was shown to fall to pre- 
vernalization levels 3 and 5 weeks after cold treatment^ ^ A gene 
loop occurring between PLC 5' and 3' regions is disrupted 
early during vernaHzation in A. thaliana. This chromosomal 
configuration change has been proposed to reveal COOLAIR 
promoter and favour its transcription^ ^ Interestingly the pre- 
vernalized PLC 5'-?>' gene loop does not reform after return to 
warm temperature. Such a stable change in conformation might 
occur at PEPl as well, and might provide a transcriptional 
memory responsible for COOLAIR priming. This could 
contribute to the increased sensitivity of older plants, which 
have already experienced a winter, to flower in response to 
inductive environmental conditions, a process often called 
competence to flower, which is particularly significant in 
perennials^^. 

Methods 

Analysis of sequences of FLC orthologues. PLC sequences from A. lyrata, 
C. rubella and T. halophila (see Supplementary Data 1 for sequences) were 
obtained using genome browsers on phytozome (www.phytozome.net) and 
checking the synteny with A. thaliana. PEPl genomic sequence was obtained from 
the BAG (accession number FJ543377.1). 

Pairwise alignments were performed using either Shuffle LAGAN^-^ and 
visualized with VISTA browser using default settings^"* or with GATAliner (using 
default settings) and visualized with GATAplotter^^. 

Plant material. A. thaliana PRI PLC was obtained by introgression of the PRI 
allele of San-Feliu-2 into Golumbia and PRI flc-3 was then recovered after EMS 
mutagenesis-^. Transgenic lines PRI flc-3 PEP lb were obtained by floral dipping^^ 
of PRI flc-3 with the Agrobacterium carrying pGAP-P£Pi^-Km plasmids (see 
Gloning) and confirmed by PGRs (primers set in Supplementary Table 3). PRI 
flc-3 PEP lb vrn2-l fca-1 plants were obtained by crossing PRI flc-3 PEP lb 18-3 
homozygote plants with vrn2-l fca-1 double mutants^^. 

A. alpina plants used in this study are from the Pajares and Dorfertal 
accessions 

Growth conditions. Flowering time was measured on 9-23 plants per genotype 
(Supplementary Table 2) grown in climate-controlled glasshouses (20 °G under 
long days, 16 h light per 8h dark). A Student's t-test was applied to determine 
whether the flowering behaviour of genotypes was statistically different from the 
control PRI^^flc3 with a P value < 0.005. Vernalization experiments were carried 
out in A. thaliana plants grown 10 days in climate-controlled glasshouses (stratified 
2 days at 4 °G), vernalized for 40 days at 4 °C under short days (8 h light per 16 h 
dark) and transferred back to climate-controlled glasshouses for 10 or 20 additional 
days. Full seedlings or rosette leaves were harvested before or after floral transition 
respectively. Fifty-six day old A. alpina plants grown in long days climate-con- 
trolled glasshouses (stratified 5 days s at 4 °C) were vernalized for up to 84 days at 
4 °C under short days and transferred back to climate-controlled glasshouses for up 
to 35 additional days in long days. To mimic two successive winter periods, some of 
the A. alpina plants vernalized as described before were transferred for another 84 
days at 4 °C under short days followed by 35 days in long days climate-controlled 
glasshouses. For each time point, one leaf of the same order was harvested from the 
main shoot on three different individuals. 

Cloning. PEPlb cloning in the pGAP-Km vector was performed using the 
recombineering method^^. The homology arms, EcoRI-Flank 1-SalI and Sall-Flank 
2-BamHI were amplified by PGR (see Supplementary Table 3 for primers) on the 
PEPl-BAC (FJ543377.1). After recombination in SW102 cells carrying PEPl-BAC, 
the pGAP-P£Pi^-Km was further extracted and transformed into Agrobacterium 
GV3101 (pPM90RK) strain. The introduced fragment was 13.4 kb comprising the 
4.4 kb PEPlb 5'-intergenic region spanning the end of PEPla upstream duplication 
to PEPlb, exon lb to 7 and a 1.8 kb downstream segment harbouring the region of 
homology at the 3' end of PEPl and PLC. 

RNA extraction. Total RNA was isolated using Trizol reagent (Invitrogen) and 
subsequently digested with DNAse I Amplification grade (Invitrogen) according to 
the manufacturer's protocols. 

3' RACE PCR. 3' RAGE experiments were performed on 3 |ig of total RNA with 
the 3' RAGE System for Rapid amplification of cDNA ends (Invitrogen) according 



to the manufacturer's protocol. PGR followed by a nested PGR were performed 
using AaGOOLAIR- specific forward primers (see Supplementary Table 3) together 
with the provided Universal Amplification Primer and the Abridge Universal 
Amplification Primer respectively. PGR products were cloned and sequenced. 

RT-qPCR and '3' tiling' RT-qPCR. Reverse transcription was performed with 
oligo (dT) primers on 1 |ig of total RNA using Superscript II Reverse Transcriptase 
(Invitrogen) according to the manufacturer's protocols. The qPGRs were per- 
formed on the Roche Light Gycler 480II instrument using the primer sets listed in 
Supplementary Table 3 and the IQ SYBR Green Supermix (Bio-Rad). Transcripts 
levels were expressed as a percentage of the accepted references genes RAN3 and 
ACTYATfor A. alpina and A. thaliana samples respectively^ ^'^^ '3' tiling' RT-qPGR 
was used to approximately determine the 3' end of AaGOOLAIR II transcripts. 
RTqPGRs utilizing the same forward primer but different reverse primers spanning 
the expected 3' end of AaGOOLAIR II (see primer sets in Supplementary Table 3) 
were performed. Primers located progressively further towards the 3' end were 
used until no amplification was detectable. PGR products were cloned and 
sequenced. 
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